Comment to "Packing Hyperspheres in High-Dimensional Euclidean Space" 
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In a recent paper [1] Skoge et al. presented new ac- 
curate numerical data on amorphous packings of hard 
spheres in space dimension d > 3. It was shown (see e.g. 
Fig. 4 in [1]) that on compressing the low density liquid 
at a constant rate 7, the pressure of the system follows 
the equilibrium pressure of the liquid up to some den- 
sity (often called glass transition density) above which 
the pressure starts to increase faster than in equilibrium, 
and diverges on approaching a value of density (pjij), 
which is called jamming density. Crystallization seems 
to be strongly suppressed by kinetic effects in dimen- 
sion d > 3 and can then be neglected in the following 
discussion. Values of 4'j{"f) have been accurately mea- 
sured in [1] as a function of 7. On the contrary, ipgij) 
is not precisely determined as long as 7 > 0: the glass 
transition is smeared and happens in a crossover region 
[^g (7)1 (7)]- However, the amplitude of the crossover 
interval seems to decrease for 7 ^ (see again Fig. 4 in 
[1], and see [2] for a recent theoretical discussion of these 
effects). 

Recently, a theory of the glass transition of hard 
spheres, that can be applied in any space dimension d, 
was developed [3, 4]. The basic idea is that the phe- 
nomenology observed at finite 7 is due to an underly- 
ing thermodynamic glass transition [5-7]: i.e. that in 
the limit 7 0, 0g (7) 4>k, the Kauzmann den- 
sity, at which an equilibrium phase transition to a glass 
phase happens. This is an idealized picture that neglects 
metastability effects due to the presence of the crystal: 
the consequences of this approximation might be impor- 
tant but cannot be discussed here, see [4] for a detailed 
discussion. However, as it seems that crystallization is 
negligible in d > 3, a comparison between the theory and 
the data of [1] is possible. In the equilibrium glass phase, 
the pressure increases faster than in the liquid phase (i.e. 
the compressibility is smaller), and diverges at a value of 
density 4)rcp = lim-y^g that we call random close 

packing (RCP) density. This is the definition of random 
close packing that can be given within our theory. How- 
ever it is affected by the metastability effects related to 
the presence of the crystal: see [4, 8] for a discussion. 
In particular, in [8] the notion of random close pack- 
ing density has been criticized and an alternative notion 
of maximally random jammed (MR J) packings has been 
proposed, see below. 

Our theory is based on standard liquid theory and on 
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0.0739 
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TABLE L Values of cpK and 4>RCp from the theory (only 
values for d < 8 are reported for brevity, values for d > 8 are in 
Fig. 1) compared with the available measured values of 4>mrj 
[1]. The last column gives the value of (f>RCP extrapolated 
from the data of [1] (see text). 



the replica trick: it takes as input the equation of state 
of the liquid phase (in practice one has to choose an ex- 
pression that describes well the liquid at low density and 
extrapolate to higher density), and gives as an output: 
i) the Kauzmann density (pK] U) the random close pack- 
ing density 4>rcp] in) the equation of state in the glass 
phase; iv) some properties of the pair correlation function 
in the glass phase, e.g. its shape close to contact; v) the 
equation of state of the metastable glass states that are 
reached for 7 > and their contribution to the entropy 
(the configurational entropy or complexity); see [3] for all 
the technical details. One of the most interesting predic- 
tions of the theory is that 2'^(j)jicp ^ dlogd for large d. 
As far as I know, this scaling has been proposed in [3] for 
the first time. The aim of this paper is to show that the 
results of [1] are fully compatible with this prediction. 

To this aim, an expression for the equation of state 
of the liquid in d > 3 is needed as input to the theory. 
The simplest choice is a generalization of the celebrated 
Carnahan-Starling equation of state to d > 3 [9]: 
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where V((/)) = gicr^) is the value of the radial distribu- 
tion function at contact, and h and are the second and 
third virial coefficients, whose exact expression is known 
[9]. The entropy of the liquid S{(t)) is obtained by inte- 
grating the exact expression 
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Given S{<j)), the random close packing density is the so- 
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FIG. 1: Plot of (fiK (open squares, obtained solving Eq. (4)), 
(jiRCP (full diamonds, Eq. (3)), and 4)mrj (full circles, nu- 
merical estimate of [1]) as a function of the dimension. Both 
(j>K and 4>RCP scale as 2''(/!> ~ dlogd for large d, while their 
distance scales as 2'*[</>i{CP — 0a'] ~ d. In the inset the same 
plot for 3 < d < 6 (compare with Fig. 6 in [1]). 
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while the Kauzmann density is the solution of 
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with Qo — 0.638... [3]. These equations can be easy 
solved numerically to get the values of (f)RCP and c/jk 
for any given value of d. The results are reported in ta- 
ble I for c? < 8, and compared with 4>mrj as reported 
in [1]. The latter quantity is the density of the maxi- 
mally random (according to some measure of "order") 
collectively jammed packings of the system, see [8] for 
the precise definition; it is estimated in [1] by the jam- 
ming density (t>j{"f) for finite but small 7 (see sec. IV in 
[1] for a detailed discussion). As 0,7(7) is expected to 
increase on decreasing 7 and 4>rcp = lini^^o </'j(7)i it 
follows that (pMRj, as estimated in [1], should be strictly 
lower than (pRCP, but close to it, consistently with the 
data in table I. A plot of (t)K and (f>RCP for d up to 19 
is reported in Fig. 1. Note that it has been suggested in 
[1] that 2'^(j)MRj ^ c\ + C2d; this scaling is not in contra- 
diction with 2'^(j>K, 2'^(j}RCP ~ dlogd because 4>mrj does 
not need to be bigger than (|)K^ even if this seems to be 
the case for d < 6. 



The very nice data for d = 4 reported in Fig. 4 of [1] 
allow for a more precise comparison of the numerical and 
theoretical results: the value of ipjij) is reported for 5 
different values of 7 = 10"^ 10""*, 10"^ 10"^, 10"'^. A 
standard procedure to extrapolate to 7 —> is to fit the 
data to a Vogel-Fulcher law: 



7(</'j) = 70IO 



-fj ^ (t)j{-f) = 4)RCP + 
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logio[7/7o] 



Such extrapolations are often ambiguous; however the 
fit is good and gives (f)RCP — 0.473, D — 0.03, 70 = 
0.45. The final result for (t)RCP differs from ~ 10% from 
the theoretical value, see table I. This is a very good 
result given the ambiguities that are present both in the 
numerical data (numerical error and extrapolation) and 
in the theory (the choice of a particular expression for the 
equation of state of the liquid that is not exact, see [10] for 
recent contributions). Note that a similar extrapolation 
is not possible in c? = 3 due to crystallization, and for 
d > 4 due to lack of numerical data. Hopefully new data 
for d > 4 will allow for a similar comparison also in this 
case. Note also that the value of (t)K ^ 0.43 we obtain 
in d = 4 seems to agree very well with the extrapolation 
of 4>g{'^) (defined roughly as the point where the curves 
leave the liquid equation of state) to 7 = in Fig. 4 of [1]. 

A more accurate comparison of the theory with the 
numerical data is possible: for instance, the theory gives 
a prediction for the glass equation of state that is close to 
the measured pressure in the glass branch for very small 
7, e.g. 7 — 10^'' in Fig. 4 of [1]; it also predicts that 
the amorphous packings are isostatic, i.e. the average 
number of contacts per sphere is z — 2d, in any dimension 
d, and this seems to be confirmed by the numerical data. 
Other properties of the packings such as the shape of the 
correlation function g(r) close to contact are predicted 
to be independent of the dimension. Unfortunately we 
are still not able to give a prediction for the splitting 
of the second peak of g{r), which seems to be strongly 
suppressed in d > 3. We hope that future work will 
address this and many other open questions [4] . The data 
reported in [1] provide a very good test of the theories 
aiming to understand the behavior of hard sphere in large 
space dimension [3, 11-13] and strongly call for further 
developments. 
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